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The dynamical fluctuations of kaon-to-pion ratios have been studied over a wide range of center- 
of-mass energies \/s. Based on changing phase space volume which apparently is the consequence 
of phase transition from hadrons to quark-gluon plasma at large ^/s, single-particle distribution 
function / is assumed to be rather modified. Varying / and phase space volume are implemented in 
the grand-canonical partition function, especially at large yi, so that hadron resonance gas model, 
when taking into account the experimental acceptance and quark phase space occupation factor, 
turns to be able to reproduce the dynamical fluctuations over the entire range of ^/s. 
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I. INTRODUCTION 



CO I Studying the collective properties of hot and dense hadronic matter is one of the main objects of heavy-ion physics. 

^H . The dependence of transport coefficients, phase structure and effective degrees of freedom on the incident energy 

and system size likely provides fruitful tools to study the collective properties. Studying fluctuations in particle 

r~| . physics of cosmic-rays dates back to early eighties of the last century [l| . The phase structure and event -by-event 

rS i' fluctuations [2|-y| have been suggested to provide comprehensive characteristics of the particle production yields in 

JL . heavy-ion collisions. They statistically variance around the averages and their scales and physically can be related to 

gj ' hadro-chemical composition of the particle source [a] and therefore to the phase transition. In this regard, examining 

(~| I the hypothesis about the equilibrium of chemical processes in the hadronic final state [g| is very much essential. The 

dynamical fluctuations of certain particle yields have been studied at SPS and RHIC energies 0-13 • Therefore, it is 

natural to study the energy-dependence of the particle yield ratios and the event -by-event fluctuations using hadron 

►^ , resonance gas model (HRG), since this model provided a good description for the thermodynamical evolution of the 

tH- ■ hadronic system below the critical temperature pi)| - ll2| and has been used to characterize the conditions deriving the 

t^ , chemical freeze-out [1^ [ij] . 

^^ ' In many aspects, the fluctuations in particle production are essential, particularly with regard to examining the 

■^ I existing statistical models [y, llOl. characterize particle equilibration [l^ and search for unambiguous signals for the 
creation of new state of matter ^^ [l3| • Within the context of the statistical models of particle production [18[ , it has 
^^ ' been concluded that none of these models is perfectly able to describe the experimental data. Using grand-canonical 
^D . ensembles and quark phase space occupation factor 7 makes the comparison relatively better. In HRG with grand- 
canonical ensembles, the dynamical fluctuations of particle yield ratios, including kaon-to-pion, have been studied 
|2(j |. To bring theoretical and experimental data very close to each others, the factor 7 played an essential role. It 
has been noticed that the fluctuations over the whole range of the available center-of-mass energies ^/s exhibit a 
. . non-monotonic behavior and 7 is varying with increasing y/s, so that SPS- are different from RHIC-data. 
j^ In present work, we recall the concept of phase space dominance suggested by Fermi |2l|, six decades ago. The 

assumption of equilibrium single-particle distribution function / and extensive thermodynamics that have been shown 
to perfectly reproduce all thermodynamic quantities, including multiplicities and fluctuations, as long as the energy 
density is not high enough to derive the hadronic system into new state of matter, quark-gluon plasma (QGP) should 
be no longer valid, when the energy density exceeds a critical value, as the case at RHIC and LHC energies. Across the 
phase transition, the symmetries and accordingly the effective degrees of freedom are likely subject of a prompt change. 
Therefore, the phase space volume, in which the microscopic states are distributed according to / turns to a subject 
of a change as weU. In present work, we implement the change of phase space volume and / in the grand-canonical 
partition function, especially at large y/s. It is not a theory describing -yi-dependence of / and/or the dynamics of 
changing phase space volume. This would be rather the role of QCD. For this purpose, Boltzmann equation integrated 
over momentum-space can alternatively be used. This has been applied on characterizing freeze-out of nucleosynthesis 
in the early universe [13] • 
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Appart from these argumentations, this idea, from the phenomenological point-of-view, does not appear from 
nowhere. While studying particle number and ratio fluctuations, there have been speculations about the possibility 
of phase space change [20| with increasing -y/s. The entropy per particle, S/N, has been analyzed over the available 
range of y/s. It was clear that S/N is able to reflect a kind of rapid modification in phase space volume. The role 
of 7, the quark phase space occupation factor, can then be carried out by S/N. It has been concluded that the 
dependence of phase space volume on energy is likely essential to access the phase transition, at large ^/s. For kaon 
production in heavy-ion collisions, especially the horn at top SPS energy, an initial partonic phase has been assumed 
beyond this energy threshold ^^. Below it, a hadronic initial state has been utilized. The dependence of kaon-to- 
pion fluctuations on both energy- and multiplicity-scaling has been analyzed, recently [23J, where the multiplicity of 
kaon- and pion-yields has been scaled out depending on the dynamical fluctuations. These studies would explain why 
kaon-to-pion ratios are thought to be interesting. Furthermore, the measured minima at top SPS energy would be 
related to an increase of entropy production and a decrease of strangeness-to-entropy ratio. In heavy-ion collisions, 
the strangeness dynamical fluctuations are of great interest, since they are sensitive to the equation of state and 
microscopic structure of strongly interacting matter created at early stage of the collisions [23| ■ 

The paper is organized as follows. The model is given in section HIl where single-particle equilibrium distributions 
of hadrons and partons arc introduced. Section IIIII is devoted to the dynamical fluctuations of kaon-to-pion ratio. 
Discussion and final conclusions are elaborated in section IIVI 



II. MODEL 

A. Single— Particle Equilibrium Distribution of Hadrons 

The grand-canonical partition function is given by Hamiltonian and baryon number operators, H and b, respectively, 
and depends on temperature T = 1//3 and chemical potential /i, 

Z{|3,V,^i) = Tr [7 exp'3(^^-")] . (1) 

It can be characterized by various but a complete set of microscopic states and therefore the physical properties of 

the quantum systems turn to be obtainable in approximation of non-correlated independent hadron resonances. The 

resonances treated as a free gas |1CHL3 . l25l . |26| are conjectured to add to the thermodynamic pressure of hadronic 

matter. This statement is valid for free, as well as, strong interactions of the resonances. It has been shown that the 

thermodynamics of str ong ly interacting system can be approximated to an ideal gas composed of hadron resonances 

with masses < 2 GeV [2^, [27[ . The main motivation of using the Hamiltonian is that it contains all relevant degrees 

of freedom of confined and strongly interacting matter. It implicitly includes the interactions that result in formation 

of new resonances. In addition, this model has been shown to provide a quite satisfactory description of the particle 

production in heavy ion collisions. 

The distribution of resonances in the available states in micro-canonical, canonical (statistical) and grand-canonical 

ensemble is subject of statistical physics. Here, we consider the occupation number basis of single-particle |28| as 

the one which is suitable to evaluate Eq. ([1]). Each macro-state is to be characterized by n g {n^} set of allowed 

TTiicrostates, their bi baryon charges and e^ = (fc^ + ?n|)^'^ energies of i-th single state of mass nii and quantized 

—> -* 2 

longitudinal conjugate momenta k. The effective mass mi ~ (r/iOf -I- fcij)^'^, where mOi is the physical mass and kt 

is the transverse momenta. In the relativistic limit kt << fc, the effective mass refers to the physical one. Obviously, 

the macro-state energy Ei is given by the sum over all possible states, i.e, Y^^niti, where the number m depends 

on whether the resonance is fermion or boson. For bosons, Ui can be very large. For fermions, Ui is subject of Pauli 

occupation principle. This summation counts all possible basis states. To make any physical treatment, one would 

need to make approximations according to the required precision. In this regard, the summation can be replaced 

by integration over the controlling variable, for instance, over e, which in turn depends on fc, where k is subject of 

quantization. Therefore, the summation over e is to be replaced by an integral over all ki in the momentum space. 

In infinite volume^ limit, V = dx^ , the phase-space integral replaces the sum over all possible 2 x n-dimensional 

phase-space volume element d'^x d^p, which apparently gives the probability of finding the particle in this infinitesimal 



^ The volume of the space volume is given by J^ Uid^xd^k, where x is the coordinate. In many cases, the geometry is cither simple or 
knowrn and, therefore, symplectic manifolds can to be assumed, wrhere a twro-form lu = cuijdx^ /2 A dx^ is fulfilled. Symplcctic forms are 
- per definition - closed and non-degenerate. 



region of phase-space, i.e, J^i — 9/i'^''^)^ Jo d^xcPp. n refers to the dimensions of the physical system. As given 
above, we used n = 3. The phase space distribution extends mainly along the large longitudinal momenta. 

In equilibrium [SOj, the average energy E = X^i'^i^i' ^^'^ particle counts N = X^i "■« of macroscopic states are 
conjectured to be know and remain constant. With respect to these two constrains, one can take the Hamiltonian 
on the system as an observable quantity. Here, we only consider the configuration which maximizes the number of 
microstates, instead of considering all configurations {ui}. The maximum count, O ~ {Ni\/Ilini\)Ilig^\ where gi 
is the degeneracy factor of the configuration rii of i-th microstatc. It is related to the global entropy S — \n{il) 
and information on the system compared to the measurements. This maximum is a formal subject to the external 
constraints on intensive variables E and N. Since the states are conjectured not to be correlated with each other, the 
chemical potential fx is vanishing. Then, the maximum number of microstates is given by solving 

Q Pi / " \ 

— (S - aN - /3E) = — i\nm + Inll'^g'^' -J^^nn^l-aN-^E] = 0, (2) 

which means that only the terms having same subscript j remain finite. The coefficients a and /? arc Lagrange multi- 
pliers in entropy maximization. Each of these multipliers basically adds some unknown amount of each independent 
constraint to the function being optimized and ensures that the constraints are satisfied. 

Ai,n,,--Ai,.^,_«_^,^_o. (3) 

Utilizing the Stirling approximation, then the occupation number. 



"J 



ffj exp(-a-/3ej), (4) 



which apparently falls off exponentially with increasing e, since, as will be shown below, 7 = exp(— a) is constant. 

The physical meaning of (3 is, as given above, the inverse of temperature [29|. Assuming two configurations n.i and 
rij, then /3 = [ln{ni/nj) + \n{gi/gj)]/{ei — ej). Therefore, /3 can alternatively be defined as the change in microscopic 
configurations per energy change. Apparently, this definition is consistent with the second law of thermodynamics. 
From thermodynamics, especially the identification of the Gibbs equation, a is another variable controlling the number 
of particles in the phase space. 

a = h.f^^lj^'] ^\nJ29^|3e^~\nN. (5) 

^ ^ i 

Then, the occupation number reads 

f N \ fNT\ 

^j = 9j l y^n ^g. ) exp(-/3ej) = .gj (^1 exp(-/3ej), (6) 

where £ = X,- giei is the summation of all possible single-particle state energies multiplied by the effective degrees 
of freedom. The multiplier a can be written as 

c^ = i^^[j^^=^^S-^^T-lnN. (7) 

It combines intensive variables, T and N with an extensive one £. The most probable state density is to be found 
by Lagrange multipliers, where one of them, a, has been expressed in term of the second one, /3, and the occupation 
numbers of the system. Apparently, a gives how the energy £ is distributed in the microstates of the equilibrium 
system and therefore, can be understood as another factor controlling the number and thus is similar to chemical 
potential at the microcanonical level. 

Besides the two constrains given previously, the conservation of baryon number ni, represents an additional constrain 
on the grand-canonical partition, Eq. ([ij. Therefore, 

Zg,iT,V,^i) = TV [exp''^-"] , (8) 

exp ( ^^ — a ] 



With these assumptions, the dynamics of the partition function can be calculated as sum over single-particle 
partition functions Z' of all hadrons and resonances. 

\nZg,{T,V,^i) = ^lnZ^,(r,V^,/i) = ^±^V^/'"fc2^fcln(l±7A,e-^'W/^) (10) 

i i ^ 

where ± stands for bosons and fermions, respectively. A^ — cxp{fii/T) is the i-th particle fugacity. 7 = exp(— a) is 
the quark phase space occupation factor. 

B. Single— Particle Equilibrium Distribution of Partons 

The Fokker-Planck equation is a well-known tool used to study the dynamics and velocity distribution of objects 
in thermal background, such as the transport properties in quark-gluon plasma |3ll 433|. Studying the stochastic 
behavior of a single object propagating with random noise, known as Langevin equation 3J, |35| . represents one way 
to solve this problem. A master equation, such as the linearized Boltzmann-Vlasov equation |36|, with the Landau 
soft-scattering approximation would give another solving method. 

As given previously, the statistical properties of an ensemble consisting of individual parton objects is given by a 
single-particle distribution function /. The probability of finding an object in infinitesimal region of the phase-space 
is directly proportional to the phase-space volume element and the distribution function describing it, where / - in 
this case - is assumed to fulfill the Boltzmann-Vlasov (BV) master equation, which is the semi-classical limit of a 
time-dependent Hartee-Fock theory through the Wigner transform of the one-body density matrix and can be used 
to study the dynamics of the constituents quarks in hadrons, 

/ + a--v,/ + fc-Vfc/ + g;-v,j = g + c. (ii) 

In rhs, the first term Q represents gain, i.e, the rate that a particle with momentum k + kt loses momentum kt due 
to reactions with the background and the second term C represents loss due to the scattering rate. Therefore, these 
two terms furnish us with details about the interaction. Details about participating partons and their reactions have 
been discussed in literature, for example in Ref. (37[. Seeking for simplicity, a two-body system has been assumed. 
Then the effective potential U has to combine the well-known Coulomb U{x -^ 0) oc 1/x and confined potentials 
U{x — )■ 00) ex 0. As assumed, the particles of interest are partons currying electric and color charges. In Ihs, standard 
position X and momentum k variables are given in the first two terms. The third term represents the dynamics of 
the charge, where k can be given by field tensor (Lorentz relativistic force) . The fourth term reflects an extension of 
phase space to include color charge, ql. 

One of the largest advantages of the transport equation, Eq. (|lip . is that the phase space distribution functions of 
partons and/or hadron resonances by test-particle distributions for these different species. In test-particle approach 
[38| , the continuous distribution is discritized to a finite number of test-particles representing the individual phase 
space cells. The test-particle itself obeys Newtonian equation-of-motion [33. Therefore, 

(12) 

(13) 
(14) 

where qc represents the chromofield exchange of quarks and gluons in the presence of the background field. F can 
be expressed, as mentioned above, in Lorentz force, qc{E + v x B) or any effecti ve p otential and higher terms can be 
replaced by the rest of the Langevin system of equations. As given in Ref. [33, 1401, BV master equation is solvable 
for 

1 ^ 
/(x,fc,i) = Y^gx{S~x,[t))gk[k~k,{t)), (15) 

where ntp = N/q^ is the total number of test-particles per quark charge. In this model, the test-particle can be the 
hadron yields, in which we are interested. The functions g^ and gk can be Heaviside or any sharply peaked distribution 
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= Fix) + ■■■, 


dqc 
dt 


= r'^'u^qbA-^ 



function [40|. In order to make the physical measurements possible, they have to be modeled [37 
Gaussian-type, 



for example as 



gx{x-Xi{t)) = W — exp{-w[x- Xi(i)]}, (16) 

where to = Xjla^ and the invariance a > 0. Positions Xi and momenta fc^ vary around x and fc, respectively. Eq. (|15|) 
can be interpreted as a measure for the deviation from equilibrium or the relaxation toward equilibrium. 

The Hamiltonian equation-of-motion of the test-particle is obtainable by inserting Eq. (|15p in Eq. pT|) . The 
function / reaches equilibrium, when the resulting probability current vanishes [33, |4l|. To illustrate this, let us recall 
the Tsallis distribution J42l | depending on the parameter q, which measures the degree of extensivity of the entropy S 
in the system. Then the single-particle distribution function reads 



/(a-,fc) 



iV 



i-:S)(i-,) 



^ 1/(1-9) 



(17) 



which obviously gives the extensive Boltzmann limit, when g — >■ 1. In Eq. (jisp . the Boltzmann limit is reachable at 
very large time i, i.e, when the number density turns to be conserved or the change in entropy S vanishes. Thus, the 
number density and distribution function would be written as 



n{x,k,t) = neq{x,k) + f{x,k,t), 
f{x,k,t) ~ feq{x,k) Q{x,k), 



(18) 
(19) 



where Q{x, k) ~ exp(2a;xfc) exponentially raises with increasing position x and momentum k of the test-particle. It 
reflects the change in the phase-space when the hadronic degrees of freedom are replaced by the partonic ones. The 
factor Q would not be too much different than q. Therefore, it can be interpreted as a measure for the non-extensivity. 
The exact relation between q and Q, especially in the context of heavy-ion collisions [43[, will be the subject of an 
upcoming work, feq {x, k) is an equilibrium normalized distribution function similar as the one assumed and utilized 

in previous subsection. In deriving Eq. ()19p . we assumed that the products of xxi and kki vanish, where Xi {ki) is 
position (momentum) of the i-th cell in the phase space volume. 

Non-equilibrium processes can be originated from different sources. Previously, we assumed that the phase transition 
is a dominant source. Because of presence of non-Markovian processes in the kinetic equation, the extreme conditions 
in ultra-relativistic heavy-ion collisions are conjectured to produce long-range color interactions and so-called memory 
effects. They likely affect the standard equilibrium distributions and even the thermalization processes toward the 
equilibrium state [4J] . We therefore suggested in Eq. ((T7| that the distribution function is to expressed out in Tsallis 
statistics. Based on distributions of partons in the microstates of the phase-space volume, a rigorous determination of 
non-extensivity is given in Eq. (|19p. Therefore, the single-particle distribution function reads in the grand-canonical 
ensemble, 



fgc{k,lJ.) = 
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1/(9-1) 



±1 



(20) 



(21) 



where in Eq. ([2T|) . the volume compensates the corresponding part of Q{x, k). This procedure gives solid reasons not 
only for the non-extensivity, but also for the assumed modification in size of system and phase space volume. Then, 
the grand-canonical partition function can be written as 



k^dk In 1 ± 



lnZg,(T,V,M) - yi±Av / k^dk In 1 ± 7 Q cxp 
i '^''^ Jo \ 

where the pressure in this grand-canonical ensemble is given as 



ei{k) - l-M 



T 
T 



-a] (q-l) 



l/(q-l)N 



T 
Pgc{T,V,i^i) = lim —\nZgc{T,V,^i). 

V— )-oo y 



(22) 
(23) 

(24) 



III. DYNAMICAL FLUCTUATIONS IN KAON-TO-PION RATIOS 



The fluctuations in particle number are given by the susceptibihty, which is the derivative of particle number in) 
wrt chemical potential /^. 
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(25) 
(26) 



When chemically relaxing system absolves the chemical freeze-out process, the hadron resonances are conjectured 
to decay either to stable particles or to other resonances. This chemical process has to be take into account in the 
particle numbers and fluctuations as follows. 



final\ 



^directx 






(27) 



(28) 



where bj^i being branching ratio for the decay of j-th resonance to i-th particle. In order to characterize the stage 
at which the chemical freeze-out takes place, we assume that the ratio s/T'^, where s is the entropy density, gets a 
constant value [ij, [ij, [2g . 

For the event -by-event fluctuations of the ratio of two particles K/n (kaon-to-pion) are [J] 



2K/r. 



((A71k)2) ^ ((An^)2) _ ^ {AriK An^) 
{ukY (riTr)^ {nii) {n-„) 



(29) 



which combine dynamical and statistical fluctuations. Third term of Eq. [^Hl counts for fluctuations from the hadron 
resonances which decay into K and tt, simultaneously. In such a mixing channel, all correlations including quantum 
statistics ones are taken into account. Obviously, this decay channel results in strong correlated particles. To extract 
statistical fluctuation, we apply Poisson scaling in mixed decay channels. Experimentally, there are various methods 
to construct statistical fluctuations [9| . Frequently used method is the one that measures particle ratios from mixing 
events. 
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Subtracting Eq. [3D] from Eg. 1291 we get dynamical fluctuations of riK/nT^ ratio. 
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(30) 



(31) 



Present model on the dynamical fluctuations and their modiflcations with varying phase space volume, especially 
when the center-of-mass energy is high enough to secure energy density able to derive the hadronic mater into QGP, 
this model is valid for any particle yield rations. The dynamical fluctuations in kaon-to-pion ratio, {K'^ + K~)/{Tr'^ + 
7r~), are given in Fig. [1] since they combine strangeness and light-boson fluctuations. They are very sensitive to 
dcconflncmcnt and chiral phase transitions, respectively. 

As introduced in previous subsections, the HRG model assumes that the hadron resonances are point-like and non- 
correlated. The dynamical fluctuations of particle yield ratios with and without Q are produced and plotted in Fig. 
[TJ It is obvious that the average multiplicity (n) , Eq. ([27| , and the dynamical fluctuations cr^ /jy^ , Eq. ([29| , are not 
strongly dependent on the volume fluctuations. Therefore, we assume - in this work - that volume fluctuations over 
the entire range of yG are minimum and thus neglected. The experimental acceptances of the different detectors have 
been taken into account. Also, the quark phase space occupation factor, 7, has been estimated. So far, we conclude 
that HRG reproduces the experimentally measure fluctuations of (K'^ + K~)/{'k~^ + tt") ratio over the entire range 
of -y/s. In generating this excellent agreement no fitting has been performed. The comparison is given in Fig. [TJ 
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FIG. 1: Dynamical fluctuations of (/i+ + A'^)/(7r+ + tt^) ratio as function of center-of-mass energy ^/s. Taking 
into account the experimental acceptance, the lower curve represents HRG results at finite value for 7. It reproduces 
perfectly the SPS data. The dashed region shows that RHIC data arc largely underestimated. The upper curve 
represents HRG results, when phase space factor Q is switched on. RHIC data perfectly matches with this curve. 
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FIG. 2: Same as in Fig. [TJ Implementing the experimental acceptances of the different detectors, keeping 7, that 
exclusively has been used to produce the curve below ^/s '-^ 17 GeV, and utilizing phase space factor Q at higher ^/s, 
result in the curve that obviously very well describes the entire experimental data set. 



IV. DISCUSSION AND CONCLUSIONS 



The scale invariant processes are conjectured to lead to intermittency. i.e, a power-law dependence of factorial 
moments F on the phase space volume, F^'^^ (SV) oc SV~^'' , where SV ~ SxSk is the phase space volume element and Sk 
is known as the intermittency exponent [45l | . The importance of phase space dimension in the intermittency approach 
has been discussed by Ochs ^^. Fluctuations and correlations in multi-particle distribution are accessible through 



intcrmittcncy. Also, information about dynamics of the deca ying: system would be available by this statistical concept 
that includes self-similarity or fractality of hadron production |47| . Also the parton shower would exhibit intermittency. 
Nevertheless, there is no guarantee that it survives until the hadronic final state |48|. In framework of Ginzburg-Landau 
theory, the intermittency in phase transition, hadrons-QGP, has been studied |49[. The intermittency evolution has 
been bound with the dynamics of QCD critical point, where it was concluded that the freeze-out profile would present 
a structure revealing traces of critical fluctuations [50| • 

Using generalization of central limit theorem applied to random cascading models. Levy stable laws have been 
introduced to describe and classify the intermittency patterns due to different kinds of phase transitions, for instance 
hadrons-QGP [5ll - l53l |. In previous work, we studied the one- and two-dimensional intermittency in Pb-Pb collisions 
at 158 AGeV/c, i.e, top SPS energy [52|- The exponents, a^, are compared with the anomalous dimensions, g, which 
has been expressed in terms of Levy stable index. It has been found that the fluctuations within a narrow rapidity 
internal fulfill the requirements of the Levy stability, referring to phase transition, i.e, positive Levy stable index. In 
the same year, an interpretation of the non-extensitivity parameter, q, in Levy and Tsallis distributions has been 
suggested |54{ . 

This discussion aims to connect the event-by-event dynamical fluctuations with the intermittency and refer to 
earlier phcnomenological results that top SPS energy has been enough to secure energy density able to dissolve the 
hadronic matter into QGP. The non-extensitivity given above in the anomalous parameter q and previously in Eq. 
(|20p is connected with departure from equilibrium and therefore has to be taken into account, when the system is 
assume to go through phase transition, as the case at top SPS energy and beyond. 

As shown in Fig. [I] the dynamical fluctuations are compared with the HRG results. The quark phase space 
occupation factor 7 is apparently not able to reproduce the whole data sent, although the experimental acceptances 
of the different detectors have been taken into considerations. In addition to the assumptions introduced in j23.[23|. 
we given a novel one. We assume that the prompt raise at ^/s ~ 17 GcV is to be understood according to modification 
in the phase space volume. To this end, we studied the single-particle distribution functions / of hadrons and partons. 
It has been concluded that / is a subject of modification, especially at large -^i. In this limit, the energy density 
available to the system turns to be high enough to cause the hadronic matter, where equilibrium / is perfectly 
able to reproduce almost all essential transport properties and thermodynamic quantities, including the dynamical 
fluctuations, to go through a phase transition to QGP. Such a phase transition apparently results in various types of 
modifications, such as symmetries and degrees of freedom. Also the configurations of microstatcs in the phase space 
volume (fixdPk and the single-particle distribution function / are not exceptions. 

We apply this model, exclusively, to the dynamical fiuctuations of (/i + + A' ^)/(7r+ + 7r^) ratios, since they combine 
both strangeness and light boson fiuctuations and have one of the largest experimental acceptances in the different 
detectors. Additionally, they are very sensitive to the deconfinement and chiral phase transitions, respectively. Fur- 
thermore, the strangeness dynamical fiuctuations are expected to survive through the mixed phase. Implementing 
7 and Q in the grand canonical partition function of HRG results in a very well description of the experimentally 
measured fiuctuations of these particle yields over the entire range of ^/s. The comparison is drawn in Fig. [2j The 
non-monotonic behavior is very well reproduced over a wide range of y/s. It is obvious, that these dynamical fluc- 
tuations non-avoidablely refer to non-extensive and non-equilibrium state of matter, that basically differs from the 
one at SPS energies. The modification of the state of matter has been combined in the factor Q{x, k) ~ exp(2a; xk) 
whose numerical value is roughly estimated and kept unchanged with changing ^/s. 
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